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SUMMARY 

This  report  describes  the  computer  program  that  was  developed  to  support  the  effort  currently 
being  conducted  at  C.  S.  Draper  Laboratory,  Inc.  to  develop  a  hybrid,  learning  augmented,  lateral- 
directional  flight  control  system.  The  output  from  the  computer  program  is  a  set  of  transfer  functions  for 
trimmed  flight  conditions  which  can  be  used  to  describe  representative  F-1 8  time  histories.  The  transfer 
functions  meet  the  Level  1  Flying  Qualities  Requirements  of  MIL-F-8785C  for  Flight  Phase  Categories  A 
and  B  of  the  F-1 8's  flight  envelope.  These  include  roll  angle,  roll  rate,  sideslip,  yaw  rate  and  lateral 
acceleration  to  lateral  stick  and  pedal  inputs.  The  computer  program  was  written  in  Visual  Basic  and  the 
coefficients  to  the  transfer  functions  are  written  to  an  output  file. 


This  report  includes  a  description  of  the  transfer  function  development,  the  computer  program 
structure  and  input/output  content  and  format.  A  sample  application  is  also  included  to  aid  the  user  in 
understanding  the  program  application. 
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1.0  INTRODUCTION 

C.  S.  Draper  Laboratory,  Inc.  is  currently  conducting  a  research  program  to  design,  develop  and 
evaluate  a  hybrid,  learning  augmented,  lateral-directional  flight  control  system  for  a  nonlinear  aircraft 
model  of  the  F-18.  This  effort  is  to  support  the  technology  program  entitled  "Aircraft  Vehicle  Technology 
Block",  Section  8.0,  "Dynamics  of  Flight",  Task  Area  8.4,  "Advanced  Flight  Controls".  The  control 
system  performance  will  be  tested  against  a  linearized  reference  model  representative  of  the  F-18  with 
Level  1  lateral-directional  flying  qualities  throughout  the  flight  envelope.  The  reference  model  will 
provide  Level  1  transient  responses  to  stick  and  pedal  inputs.  This  report  describes  the  approach,  basic 
assumptions,  and  the  computer  code  used  to  develop  the  reference  model. 


2.0  APPROACH 

The  objective  of  this  effort  was  to  develop  a  lateral-directional  reference  model  that  would  provide  time 
histories  that  demonstrate  Level  1  lateral-directional  flying  qualities  throughout  the  flight  envelope  of  the 
F-18  modeled  at  C.  S.  Draper  Laboratory,  Inc.  The  approach  taken  was  to  provide  simplified  transfer 
functions  through  pole-zero  placement  that  are  representative  of  the  F-18.  The  simplified  transfer 
functions  were  designed  to  meet  the  Level  1  lateral-directional  flying  qualities  requirements  outlined  in 
MIL-F-8785C  (Reference  1)  for  Flight  Phase  Categories  A  and  B.  Only  those  requirements  which  could 
be  represented  by  transfer  functions  were  addressed  in  this  effort.  The  reference  transfer  functions  were 
then  used  to  generate  the  required  time  histories  to  different  stick  and  pedal  inputs  at  any  trimmed  flight 
condition.  For  simplification  it  was  assumed  that  the  reference  model  could  be  represented  by  the  classical 
lateral-directional  3  degree-of-ffeedom  (DOF)  equations  of  motion  (EOM),  as  described  in  Reference  2. 
With  this  assumption  there  was  no  coupling  between  the  longitudinal  EOM  and  the  lateral-directional 
EOM.  By  using  the  classical  lateral-directional  3  DOF  EOM,  the  modal  characteristics  and  the  resulting 
transfer  functions  were  described  purely  by  the  roll  mode  time  constant  (tr),  the  spiral  mode  time  constant 
(Ts);  and  the  Dutch  roll  frequency  (coj)  and  damping  (q^). 

The  poles  or  the  modal  characteristics  of  the  transfer  functions  were  taken  from  Reference  3  as  a 
function  of  Mach  and  altitude.  In  those  instances  where  the  modal  characteristics  do  not  meet  the  Level  1 
Flying  Qualities  Requirements  of  MIL-F-8785C,  the  values  were  modified  to  the  appropriate  maximum  or 
minimum  values.  This  was  based  on  the  assumption  that  the  new  control  system  would  augment  the  F-18 
model  dynamics  to  meet  MIL-F-8785C  requirements.  For  simplification  the  spiral  mode  time  constant  was 
set  to  zero  throughout  the  flight  envelope.  Although  most  of  MIL-F-8785C  requirements  could  be  satisfied 
through  the  roll  response  to  lateral  stick  force  transfer  function,  the  intent  of  this  effort  was  to  provide 
multiple  Level  1  reference  time  histories  for  the  hybrid  lateral-directional  flight  control  system.  The 
transfer  functions  gains  and  zeros  were  either  placed  or  calculated  to  provide  Level  1  Flying  Qualities.  The 
values  for  the  numerator  gains  and  zeros  depended  on  the  transfer  function  being  evaluated.  The  following 
paragraphs  will  describe  how  these  transfer  functions  numerators  and  gains  were  computed. 
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2.1  ASSUMPTIONS 


This  section  presents  the  basic  assumptions  used  to  simplify  the  equations  of  motion  in  order  to  generate 
the  transfer  functions  calculated  by  the  computer  program  (F18LDTF).  These  are  standard  assumptions  as 
applied  in  Reference  2. 


Assumption  1. 
Assumption  2. 
Assumption  3. 
Assumption  4. 
Assumption  5. 

Assumption  6. 


Assumption  7. 
Assumption  8. 

Assumption  9. 


Assumption  10 
Assumption  11. 


The  airplane  is  assumed  to  be  a  rigid  body. 

The  earth  is  assumed  to  be  an  inertial  frame  of  reference. 

The  mass  and  mass  distribution  of  the  airplane  are  assumed  to  be  constant. 

The  XZ  plane  is  assumed  to  be  the  plane  of  symmetry. 

The  disturbances  from  the  steady  flight  conditions  are  assumed  to  be  small 
enough  that  the  small  angle  approximations  are  applicable. 

The  steady  lateral  trim  conditions  are  assumed  to  be  such  that  initial  roll  (p)  and 
yaw  (r)  rates,  roll  angle  (tp)  and  sideslip  angle  (  P)  are  zero  and  the  longitudinal 
forces  and  moments  due  to  lateral  perturbations  about  such  trim  conditions  are 
assumed  negligible. 

The  air  flow  is  assumed  to  be  quasi  steady. 

The  variation  of  atmospheric  properties  such  as  density  are  assumed  to  be 
negligible. 

The  earth  is  assumed  to  be  flat  such  that  the  effects  associated  with  the  rotation 
of  the  vertical  axis  of  the  airplane  while  moving  in  the  inertial  frame  are 
assumed  negligible.  In  addition  the  airplane  is  trimmed  for  straight  and  level 
flight  with  pitch  rate  (q)  =  0 

The  equations  of  motion  are  based  on  a  stability  axis  system. 

In  the  steady  flight  condition,  the  flight  path  of  the  airplane  is  assumed  to  be 
horizontal. 


In  addition  the  following  auxiliary  relationships  were  applied; 


p=  cp 

(1) 

r=  vj/ 

(2) 

aycg  =UoP-g‘P  +  Uor 

(3) 

2.2  DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

The  assumptions  and  auxiliary  equations  cited  in  the  previous  section  permits  the  decouplmg  of  the 
longitudinal  and  lateral-directional  EOM.  In  addition  the  3-DOF  Dutch  roll  approximation  was  applied 
since  it  complemented  the  assumption  that  Tj  =  0  throughout  the  flight  envelope.  The  assumptions  used  in 
the  3-DOF  Dutch  Roll  approximations  are: 

Assumption  12.  The  rolling  acceleration  due  to  yaw  rate,  L^r,  is  negligible. 

Assumption  13.  The  yawing  acceleration  due  to  roll  rate,  NpP,  is  negligible. 

g(p 

Assumption  14.  The  contribution  of  gravity  term  to  the  side  force  equation,  is  negligible. 

'-’o 

The  resulting  equations  for  the  F-18  modeled  in  the  computer  program  are  as  follows: 
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(s-Y,)p(s)  +  r(s)  =  (Ys*^)6a(s)  +  (Y5*p5,(s)  +  (Ys*^)5,(s)  (4) 

-L^P(s)  +  s(s-Lp(p(s)  =  (L^3)6a(s)  +  (4^)53{s)  +  (L^^)8r(s)  (5) 

-N^p(s)  +  (s-N')r{s)  =  (N^^)5s{s)  +  (N^p5s(s)  +  (N^^)5,(s)  (6) 

The  transfer  functions  were  obtained  by  solving  the  Laplace  transformed  perturbation  equation  for  the 
output  variable  of  interest  with  all  other  inputs  considered  to  be  zero.  The  equations  for  the  lateral- 
directional  responses  to  controller  inputs  were  defined  as  follows: 


P(s)/5(s)  =  (Ap  s  -1-  Bp)/Alat-dir 

(7) 

r(s)/5(s)  =  (Aj.  s  +  Br)/Alat-dir 

(8) 

(p(s)/5(s)  =  (Acp  s^  +  B(ps  +  Ccp)/Alat-dir 

or  =  A(p(s2  2q^  «o<pS  +  ffi<p2)/Alat-dir 

(9) 

p(s)/5(s)  =  s(A(p  s^  +  Bcps  +  C(p)/Alat-dir 

(10) 

Alat-dir  =  s(As^  +  Bs-  +  Cs  -t-D) 

or  =  s(s  +  1/t|^)(s2  +  2qd®dS  +  cod^) 

(11) 

where: 


A,  B,  C,  D  =  coefficients  of  the  transfer  functions 
Alat-dir  =  denominator  of  the  transfer  function  . 

The  roll  angle,  roll  rate  and  sideslip  response  to  aileron  command  (cp/Sg,  p/Sg,  and  p/Sg  respectively) 
were  defmed  by  pole-zero  placement  to  meet  the  Level  1  lateral-directional  flying  qualities  requirements  of 
MIL-F-8785C.  This  was  done  by  using  F-18  data  (Reference  3)  to  define  the  denominator  of  the  transfer 
function  and  set  the  numerator  zeros  and  gains  to  satisfy  MIL-F-8785C.  The  remaining  transfer  functions 
were  defined  from  F- 1 8  aerodynamic  data  from  Reference  4.  Table  1  provides  the  coefficients  for  the 
numerators  of  the  transfer  functions  in  terms  of  the  primed  dimensional  stability  derivatives  as  defined  in 
Reference  2.  In  addition  many  of  MIL-F-8785C  requirements  are  stated  in  terms  of  responses  to  lateral 
stick  force  or  pedal  force  inputs.  In  order  to  model  the  transfer  functions  as  responses  to  force  inputs, 
simplified  linear  relations  between  force  inputs  and  control  deflections  were  derived  and  applied  as  a  gain 
to  the  transfer  functions.  By  using  this  method  the  roll  and  yaw  contributions  due  to  the  stabilator  were 
combined  with  the  aileron  inputs  to  yield  the  yaw  rate  response  to  lateral  stick  force  inputs. 


TABLE  1 

LATERAL-DIRECTIONAL  NUMERATOR  TRANSFER  FUNCTION  COEFFICIENTS 


1 

A 

B 

c 

Ys, 

Y6,(-N;-(Ni,/Y5\)) 

m 

Ns, 

N6,(-Y,+(Ys;/Ni,)Np) 

Or 

La, 

4,(-Yv-n;  +  (Ys\/4,)4) 

Na, 

Ni,(-Y,+.(Y5\/NJ,)N[,) 

mm 

Na, 

N6,(-Y,+(Ys;/N^^)N^) 
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2.3  MODAL  CHARACTERISTICS 

The  modal  characteristics  are  defined  as  the  roots  of  the  characteristic  equation  of  the  transfer  function 
denominator.  Equation  (11)  describes  the  modal  characteristics  of  the  reference  model  used  in  this 
computer  code  and  was  used  to  satisfy  the  following  MIL-F-8785C  requirements: 

3 .3 . 1 . 1  Lateral-Directional  Oscillations  (Dutch  Roll) 

3. 3. 1.2  Roll  Mode 

3 .3 . 1 .3  Spiral  Stability 

As  stated  earlier  the  spiral  mode  time  constant  (xj  was  set  to  zero  throughout  the  flight  envelope  which 
meets  the  flying  qualities  specification  requirement  in  MIL-F-8785C  Paragraph  3.3. 1.3  since  Xg  is  neutrally 
stable.  The  values  for  the  remaining  roots  of  the  characteristic  equation  were  from  Reference  3  as  a 
function  of  Mach  and  altitude  to  provide  values  representative  of  the  F-18.  The  Dutch  roll  frequency  and 
damping  were  tested  against  the  specification  minimum  values  of  Oj  >  1.0  rad/sec  and  (;(i>0.4  (MIL-F- 
8785C  Paragraph  3.3. 1.1).  Iftoj  and  did  not  meet  the  minimum  values  then  the  program  set  the  data 
to  those  minimum  values.  If  Xr  exceeded  the  maximum  specification  value  of  1 .0  seconds  (MIL-F- 
8785C  Paragraph  3.3. 1.2)  the  program  sets  Xr=  1.0.  Although  no  minimum  value  is  called  out  in  the 
specifications,  past  experience  (Reference  5)  has  shown  that  designers  should  avoid  too  low  an  augmented 
value  for  xr  to  prevent  roll  ratcheting.  Since  the  F-18  by  design  has  relatively  low  values  for  xr,  a 
minimum  of  0.3  seconds  was  selected  to  prevent  potential  roll  ratcheting  in  the  hybrid  control  system. 


2.4  ROLL  REQUIREMENTS 

The  roll  angle  and  roll  rate  response  to  lateral  stick  force  inputs  were  used  to  satisfy  the  following 
MIL-F-8785C  requirements: 

3. 3.2.2  Roll  Rate  Oscillations 

3. 3.2.2. 1  Additional  Roll  Rate  Requirement  For  Small  Inputs 
3. 3. 2. 3.  Bank  Angle  Oscillations 

3.3.4. 1  Roll  Performance  for  Class  IV  Airplanes 

3. 3.4. 1.3  Roll  Response 

3. 3.4.4  Linearity  of  Roll  Response 

The  roll  response  to  lateral  stick  force  inputs  (cp  /  Fj)  and  the  roll  rate  response  to  lateral  stick  force 
inputs  (p  /  Fj)  were  represented  by  the  following  equations: 


cp  _  +2^cp“(pS  +  “^) 

Fs  (s)(s-hl/xR)(s^+2gG)dS  +  (Od) 

p  _  Lpa(s)(s^ -i-2Qp(0(pS  +  (0^) 

Fs  (s)(s+l/xR)(s^+2gQ)dS  +  a)d) 


(12) 


(13) 


where: 


Lpa  =  command  gain  to  stick  force 

To  meet  the  roll  oscillation  requirement  stated  in  paragraphs  3.3.2.2,  3.3.2.2.1,  and  3.3 .2.3  of  MIL-F- 
8785C  the  complex  roots  of  the  numerator  and  denominator  were  canceled  (coj  =  co^,  and  =  q^).  By 
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canceling  the  complex  roots  the  Dutch  roll  mode  was  not  excited  and  thus  inherently  satisfied  the  above 
requirements. 

The  resulting  equation  for  the  roll  rate  response  was: 


P  L'pa 
Fs  (s  +  1/tr) 

The  term  L'p^  was  then  calculated  from  the  steady-state  roll  equation  (p^g)  where: 

Pss  =LFa(TR)Fs 

By  rearranging  Equation  (15)  the  computer  program  solves  for  L'p^  using  Equation  16. 


(14) 


(15) 


Pss 

TrFs 


(16) 


The  data  for  p^j  was  taken  from  Reference  6,  which  contained  data  for  predicted  maximum  steady  state 
roll  as  a  function  of  Mach  and  altitude.  The  lateral  stick  force  input  was  modeled  to  reflect  that  for  Class  IV 
aircraft  aileron  eflectiveness  which  is  nonlinear  with  deflection  (i.e.  large  deflections  produce  an 
incrementally  larger  rolling  moment  than  do  small  deflections).  This  was  achieved  by  modeling  a 
nonlinear  stick  shaping  curve  within  the  region  recommended  in  Reference  5  and  is  shown  in  Figure  1. 
Reference  5  indicated  that  designing  within  this  region,  excessive  roll  sensitivity  can  be  avoided  while 
maintaining  adequate  roll  control  power. 


Figure  1.  Design  Roll  Command  Shaping 


5 


NAWCADWAR-95030-4.3 


Paragraph  3.3.4. 1  is  a  measure  of  time  to  change  bank  angle  (cp)  within  a  specified  time  as  shown  in  Table 
2.  Table  3  contains  the  speed  range  definitions  for  Level  1  Class  IV  airplanes.  Paragraph  3.3.4. 1.3  states 
that  "Stick-controlled  Class  IV  airplanes  in  Category  A  Flight  Phase  shall  have  a  roll  response  to  roll 
control  force  not  greater  than  ISdegreesin  1  second  per  pound  for  Level  1...".  To  meet  paragraphs  3.3.4.1 
and  3.3.4. 1 .3  of  MIL-F-8785C  Lpa  is  substituted  into  equation  (12)  since  these  two  requirements  are  a 
function  of  roll  or  bank  angle.  The  computer  program  calculates  the  roll  response  for  the  given  time  span 
and  force  input  and  tests  the  resulting  values  against  the  specification  requirements  for  stick  sensitivity  and 
roll  performance.  If  the  value  for  Lpa  did  not  satisfy  those  two  requirements  the  value  for  Lpa  was  either 
increased  or  decreased  so  that  the  specification  requirements  were  satisfied. 


TABLE  2.  LEVEL  1  ROLL  PERFORMANCE  FOR  CLASS  IV  AIRPLANES 
FLIGHT  PHASE  CATEGORY  A 

TIME  TO  ACHIEVE  THE  FOLLOWING 
BANK  ANGLE  CHANGE  (sec) 


SPEED  RANGE 

30° 

O 

O 

kO 

90° 

VL 

1.1 

L 

1.1 

M 

1.3 

H 

1.1 

TABLE  3.  LEVEL  1  AIRSPEED  RANGE  DEFINITIONS  FOR 
CLASS  IV  AIRPLANES 


SPEED  RANGE 
SYMBOL 

EQUIVALENT  AIRSPEED  RANGE 

VL 

Vo„,in^V<V^j„+20KTS 

L 

V^in+20KTS<’><V<1.4Vn,i„ 

M 

■4V„,in<V<0.7V^ax^2) 

H 

(1)  or  Vg^.^  whichever  is  greater 

(2)  or  V„  whichever  is  less 
^  '  Omax 


The  requirements  of  3.3.4.4,  Linearity  of  Roll  Response  would  inherently  be  met  by  the  low  order 
representation  of  roll  and  roll  rate  transfer  functions. 


2.5  SIDESLIP  REQUIREMENTS 

The  sideslip  (  p)  response  to  lateral  stick  force  input  and  rudder  pedal  force  input  was  used  to  satisfy 
the  following  MIL-F-8785C  requirements: 

3. 3. 2.4  Sideslip  Excursions 

3. 3.2.4. 1  Additional  Sideslip  Requirement  for  Small  Inputs 

3. 3.2. 5  Control  of  Sideslip  in  Rolls 

3.3.6. 1  Yawing  Moments  in  Steady  Sideslips 
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Paragraphs  3.3.2.4,  3.3.2.4.1  and  3.3.2.5  were  satisfied  by  assuming  that  the  control  system  is  capable 
of  eliminating  any  sideslip  generated  by  the  lateral  stick  and  basically  setting  the  sideslip  response  to  zero 
(0) .  An  example  procedure  is  outlined  in  Reference  5  that  designs  a  stick  crossfeed  to  the  rudder  pedals 
which  would  ideally  eliminate  the  sideslip  generated  by  the  stick  input. 

Paragraph  3.3.6. 1  was  satisfied  through  the  sideslip  response  to  pedal  force  input  (P/Fp)  transfer 
function.  In  this  effort  the  gain  on  the  pedal  to  rudder  deflection  was  set  so  that  the  right  pedal  (or  positive 
input)  would  generate  a  left  (or  negative)  rudder  deflection  and  vice-versa. 

The  remaining  transfer  functions:  roll  and  roll  rate  to  pedal  force  input  (tp  /  Fp  and  p  /  Fp)  ,  yaw  rate  to 

stick  and  pedal  force  inputs  (r  /  Fj  and  r  /  Fp)  and  lateral  acceleration  at  the  center  of  gravity  to  stick  and 

pedal  force  inputs  (ay^g  !  Fj  and  ay^^  /  Fp)  were  included  as  a  matter  of  completeness  but  were  not  tailored 

to  meet  any  specification  requirements.  The  coefficients  of  the  transfer  functions  were  calculated  as  shown 
in  Table  i .  The  equation  for  lateral  acceleration  at  the  center  of  gravity  to  stick  and  pedal  inputs  were 
computed  according  to  Equation  (3).  Table  4  lists  the  equations  for  the  dimensional  stability  derivatives 
and  the  primed  derivatives  used  in  Table  1.  The  non-dimensional  stability  derivatives  are  provided  as 
nonlinear  table  data  as  a  function  of  Mach  and  angle  of  attack  in  Appendix  C. 

TABLE  4.  LATERAL-DIRECTIONAL  STABILITY  DERIVATIVES 


Quantity 

Definition 

Units 

Yv 

psu 

2m  yp 

1 

sec 

Ye 

psu^ 

2m  ys 

ft 

sec^ 

Y* 

psu_ 

2m  ys 

1 

sec 

Np 

pSU^b 

2Izz  "P 

1 

sec^ 

Nr 

pSUb^^ 

1 

41 

^*ZZ 

sec 

Ng 

psu^b 

21^  "P 

1 

sec^ 

N| 

Ni+O'^JLi 

l-(4/IxxIzz) 

Lp 

pSUb^ 

21  'P 

1 

sec^ 

Lr 

psub^ 

41 

1 

sec 

Le 

psu^b 

2Ixx  'P 

1 

2 

sec 

Li 

Li+(‘’%x)Ni 

l-dL/IxxIzz) 
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3.0  PROGRAM  STRUCTURE 

The  computer  program  F18LDTF.EXE  is  a  menu  driven  program  written  in  Visual  Basic  Version  2.0 
and  runs  in  the  Microsoft  Windows'^'^  Environment.  The  program  in  broken  into  two  main  sections, 
F18LDTF.BAS  and  F18LDTF.FRM  shown  in  Figure  2.  The  following  bold  face  sections  describe  each 
main  section  as  well  as  code  and  subroutines  corresponding  to  each  menu  option  in  alphabetical  order.  A 
program  listing  is  in  Appendix  A.  Modeled  aircraft  specifications  and  the  permissible  flight  envelope  are 
in  Appendix  B.  The  plots  for  the  nonlinear  data  are  shown  in  Appendbc  C.  In  addition  to  the  program  and 
data  table,  an  input  file  with  trim  Mach,  altitude  and  angle-of-attack  is  required  to  run  the  program.  A 
sample  input  file  is  in  Appendix  D.  Note  that  it  includes  a  case  outside  of  the  permissible  flight  envelope 
to  test  for  potential  input  errors.  The  output  file  lists  the  coefficients  of  the  transfer  functions  for  roll  rate, 
roll  angle,  sideslip,  yaw  rate  and  lateral  acceleration  to  lateral  stick  and  pedal  inputs  for  each  flight 
condition  evaluated.  A  sample  output  file  corresponding  to  the  input  file  of  Appendix  D  is  in  Appendix  E. 

calc_rho  -  Subroutine  that  calculates  the  air  density  as  a  function  of  altitude.  The  equations  are  based 
on  standard  atmospheric  relationships. 

Files  -  Menu  caption  "Files".  This  subroutine  activates  the  sub-menu  options  open_tig_Click, 
output_file_Click  and  flight_cond_Click.  The  intent  of  this  option  is  to  allow  the  user  to  use  and  define 
different  input  and  output  data  as  appropriate. 

flight_cond_Click  -  Menu  caption  "Open  Flight  Conditions".  Opens  a  user  defined  input  file. 

Operates  the  same  as  a  Windows™  file  open  command.  The  input  file  format  is  that  the  first  line  contains 
the  run  title.  The  following  lines  lists  trim  Mach,  altitude  in  feet,  and  angle-of-attack  in  radians  separated 
by  commas  for  each  flight  condition.  There  is  no  limit  to  the  number  of  flight  conditions  to  be  listed. 

Form_Load  -  Loads  initial  data  and  sets  default  values  at  the  initial  execution  of  the  program.  Also 
sets  the  error  output  file  for  the  Dynamic  Link  Library  (DLL)  VBREAD  to  "c:\fl  81dtf\error.out". 

F18LDTF.BAS  -  Declares  the  global  variables  and  arrays,  constant  values  and  initializes  the  DLL 
TOPEN  and  TLOOK.  TOPEN  and  TLOOK  were  modified  from  Reference  7  to  provide  a  look-up  routine 
for  the  nonlinear  table  data. 

F18LDTF.FRM  -  Loads  initial  data  and  sets  default  values  at  the  initial  execution  of  the  program. 

This  is  the  code  which  also  acts  as  the  driver  to  execute  the  various  menu  driven  subroutines  to  read  the 
data  tables,  input  files,  calculate  the  transfer  functions  and  write  the  transfer  functions  to  an  output  file. 
Each  menu  option  is  activated  by  placing  the  mouse  pointer  over  the  menu  option  and  clicking  the  left 
mouse  button.  Figure  3  shows  the  initial  display  when  the  program  is  activated. 

Exit_Click  -  Menu  caption  "Exit".  Exit  program. 

mach_to_vel  -  Subroutine  that  calculates  velocity  in  ft/sec  as  a  function  of  Mach  and  altitude.  The 
equations  are  based  on  standard  atmospheric  relationships. 

open_tlg_Click  -  Menu  caption  "Open  Tigs  File".  This  option  operates  the  same  as  a  Windows™  file 
open  command  with  the  default  nonlinear  data  table  "c:\fl81dtf\fl81d.dat".  The  data  table  is  to  be  used 
later  in  the  Dynamic  Link  Library  (DLL)  VBREAD. 

output_file_Click  -  Menu  caption  "Output  File  Name".  This  option  operates  the  same  as  a  Windows™ 
file  open  command  with  the  default  output  file  name  "c:\fl81dtf\fl  8tf  out".  The  user  may  change  the  name 
as  appropriate. 

print_t£)c  -  Subroutine  that  prints  the  Laplace  transfer  function  in  descending  order  of  Laplace  operator 
(s)  to  the  output  file.  The  x  represents  the  order  of  the  transfer  function. 
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Run_Click 

Generate  Transfer  Functions 


Figure  3.  F1 8LDTF  Initial  Screen 
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run_Click  -  Menu  caption  "Generate  Transfer  Functions".  This  code  does  the  actual  computation  of 
the  transfer  functions  for  each  flight  condition  in  the  input  file  and  writes  the  results  to  the  output  file 
named  in  output_file_Click.  The  default  output  file  name  is  "c:\fl  81dtf\fl  Stf.out". 

set_comma_Click  -  Menu  caption  "Comma  Delimiters".  This  subroutine  sets  the  delimiters  of  the 
coefficients  for  the  transfer  functions  in  the  output  file  to  a  comma  between  values. 

set_space_Click  -  Menu  caption  "Space  Delimiters".  This  subroutine  sets  the  delimiters  of  the 
coefficients  for  the  transfer  functions  in  the  output  file  to  a  single  space  between  values.  This  is  the  default 
setting. 

set  tab  Click  -  Menu  caption  "Tab  Delimiters".  This  subroutine  sets  the  delimiters  of  the  coefficients 
for  the  transfer  functions  in  the  output  file  to  tab  positions  0, 20, 40,  60,  and  80  as  appropriate  to  the  order 
of  the  transfer  function. 

Set_Up  -  Menu  caption  "Set-up".  This  subroutine  activates  the  sub-menu  options  set  space, 
set  comma,  and  set  tab.  The  program  defaults  to  a  space  delimiter.  This  option  allows  the  user  to 
customize  the  output  file  by  defining  the  delimiters  to  be  used  in  printing  out  the  transfer  functions. 

shell_Click  -  Menu  caption  "Shell".  This  subroutine  allows  the  user  to  exit  the  program  to  run  other 
executable  files  and  then  return  to  the  program. 

speed_cat  -  Subroutine  that  determines  the  speed  category  (VL,  L,  M,  and  H)  in  which  to  test  MIL-F- 
8785C  Requirement  3.3.4. 1 .  Subroutine  inputs  are  Mach,  altitude  and  velocity.  The  subroutine  calls  vmin 
and  vmax  to  determine  the  speed  limits  for  the  given  altitude  as  defined  in  the  flight  envelope  in  reference 
6.  These  limits  are  then  used  to  determine  the  speed  category  of  the  trim  velocity. 

stick_gain  -  Subroutine  that  calculates  the  linearized  gains  between  the  lateral  stick  force  input  and  the 
aileron  and  stabilator  deflections.  The  gains  are  based  on  the  values  from  References  3  and  4. 

test_drdamp  -  Subroutine  that  tests  nonlinear  table  value  for  Dutch  roll  damping  (  q^).  If  qj  <  0.4 
then  the  code  sets  q^  =  0.4. 

test_envelope  -  Subroutine  that  tests  whether  or  not  the  trim  conditions  are  within  the  permissible 
flight  envelope  of  the  F-18  modeled  in  this  program.  The  program  tests  against  the  envelope  as  a  function 
of  Mach  and  altitude.  It  does  not  test  for  stall  angle-of-attack. 

test_freq  -  Subroutine  that  tests  nonlinear  table  value  for  Dutch  roll  frequency  ( (O  ^ ).  If  (o  j  <  1 .0 
rad/sec  then  the  code  sets  (o  (j  =  1.0  rad/sec. 

test_phi  -  Subroutine  that  tests  the  roll  sensitivity  and  time  to  roll  requirements.  If  the  time  to  roll 
requirement  is  not  met  then  the  numerator  gain  is  increased  by  1 0%  until  the  requirements  is  met.  If  the 
roll  sensitivity  is  too  high  the  numerator  gain  is  set  to  the  maximum  requirement  value. 

test_stick  -  Subroutine  that  determines  the  maximum  stick  force  input  to  as  a  function  of  pgj.  The  stick 
force  is  the  same  as  the  design  curve  shown  in  Figure  1 . 

vmax  -  Subroutine  that  computes  the  maximum  speed  as  a  function  of  altitude  for  the  permissible 
flight  envelope  for  the  F-1 8  at  a  gross  weight  of  42097  lbs. 

vmin  -  Subroutine  that  computes  the  minimum  speed  as  a  function  of  altitude  for  the  permissible  flight 
envelope  for  the  F-18  at  a  gross  weight  of  42097  lbs. 
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4.0  PROGRAM  OPERATION 


4.1  SYSTEM  REQUIREMENTS 


The  program  was  designed  in  the  Microsoft®  Visual  Basic  ™  Version  2.0  programming  system  for 
Windows™  Ver.  2.0  higher.  The  program  runs  from  mouse  inputs  and  requires  1 .2  MB  disk  space  to  load 
the  distribution  disk  onto  the  host  system.  The  distribution  disk  contains  the  executable  code,  the  program 
source  code,  the  custom  controls  (denoted  with  the  extension  "VBX")  required  to  run  the  program,  the 
dynamic-link  library  (DLL)  executable  and  source  codes  for  VBREAD  that  executes  the  table  look-up 
routine  for  the  F-18  data  tables  used  in  the  program  and  a  sample  input  file  (FlSTRIM.fNP)  and  output  file 
(F18TF.OUT). 


4.2  INSTALLATION 

The  Setup  program  will  install  the  program  and  source  code  in  a  directory  called  "fl81dtf'.  For 
reference  purpose,  it  will  be  assumed  that  the  floppy  drive  designation  is  A,  the  hard  drive  designation  is  C 
and  the  user  is  in  the  Windows’’’^  environment.  To  install  the  program  insert  the  floppy  disk  into  the  A 
drive.  In  the  Program  Manager,  choose  Run  from  the  File  Menu,  and  then  type: 

A: setup  <RETURN> 

The  Setup  program  will  create  and  copy  the  following  files  into  the  C:\F1 8LDTF\  directory: 

F18LDTF.BAS 

F18LD.DAT 

F18LDTF.EXE 

F18LDTF.FRM 

F18LDTF.MAK 

F18TF.OUT 

F18TRIM.INP 

ERROR.OUT 

VBREAD.DLL 

VBREAD.FOR 

VBREAD.DEF 

Into  C:\WINDOWS\SYSTEM\  directory  the  following  custom  control  files  will  be  copied: 

ANIBUTON.VBX 

CMDIALOG.VBX 

GAUGE.VBX 

GRAPH.VBX 

GRID.VBX 

KEYSTAT.VBX 

MSCOMM.VBX 

MSMASKED.VBX 

OLECLIEN.VBX 

PICCLIP.VBX 

SPIN.VBX 

THREED.VBX 

VBRUN200.DLL 
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4.3  PROGRAM  EXECUTION 

The  following  instructions  will  describe  how  to  invoke  "FISLDTF"  and  run  a  sample  session.  To 
invoke  the  executahle  program  go  to  the  Program  Manager  in  the  Windows  environment.  Choose  Run 
from  the  File  Menu  and  type: 

C:\F18LDTFMF18LDTF  <RETURN> 

Figure  3  shows  what  will  appear  on  the  screen  at  program  execution.  To  run  the  program  use  the  mouse  to 
click  on  the  various  menu  options.  The  following  instructions  will  describe  the  recommended  order  in 
which  to  exercise  the  various  menu  options. 

1.  Move  the  mouse  to  "Files"  and  select  the  option  "Open  Tigs  File".  The  default  file 
"c:\fl81dtf\fl81d.dat"  will  appear.  Click  the  "OK"  button.  When  the  data  file  is  read  in,  the 
message  "DLL  LIB  TOPEN  COMPLETE"  will  appear  on  the  screen  in  the  "STATUS"  text  box  . 

2.  Move  the  mouse  to  "Files"  and  the  option  "Open  Flight  Conditions"  and  click  on.  This  option 
operates  the  same  as  a  Windows''^^  file  open  command.  The  sample  input  file  is  in 

"c:\fl  81dtf\fl  8trim.inp".  The  user  may  create  their  own  input  file  for  transfer  function  file 
generation.  Click  the  "OK"  button  when  the  appropriate  input  file  is  highlighted.  In  the 
"STATUS"  text  box  a  message  will  appear  stating  input  flight  condition  file  is  open. 

3.  Move  the  mouse  to  "Set  Up".  The  default  delimiter  between  the  coefficients  of  the  transfer 
function  is  a  space.  If  the  user  requires  a  different  delimiter  the  user  may  select  the  option  for  a 
comma,  space  or  tab  delimiter.  A  message  will  appear  in  the  "STATUS"  text  box  stating  which 
delimiter  was  selected.  If  the  space  is  a  satisfactory  delimiter  then  the  user  may  skip  this  option. 

4.  Select  the  "Generate  Transfer  Functions"  option.  While  the  program  is  reading  in  the  flight 
conditions  and  generating  the  transfer  functions  the  message  "Generating  Transfer  Functions"  will 
appear  in  the  "STATUS"  text  box.  In  the  box  below  the  "STATUS"  text  box  the  message  "Case 
Number  x"  will  appear  for  each  flight  condition  case.  Upon  completion,  the  message  "Generation 
of  Transfer  Functions  Complete"  and  "Transfer  Functions  Saved  to  File  FI  8TF.OUT"  will  appear 
in  the  "STATUS"  text  boxes. 

5.  Select  "Exit"  to  leave  the  program. 

6.  If  the  user  desires  to  run  additional  executable,  batch  or  command  files,  the  "Shell"  option  may 
be  invoked.  This  is  useful  if  the  user  would  like  to  view  or  edit  the  output  file  prior  to  exiting  the 
program. 
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5.0  CONCLUSIONS 

This  report  documents  the  development  of  a  computer  program  that  generates  Level  1  lateral-directional 
transfer  functions  for  a  Class  IV  airplane  that  meet  requirements  of  MIL-F-8785C  given  a  trim  Mach, 
altitude  and  angle-of-attack  flight  conditions.  The  data  used  to  generate  the  transfer  functions  were 
representative  of  an  F- 18  and  will  support  the  effort  being  done  at  C.  S.  Draper  Laboratory,  Inc.  to  develop 
a  hybrid,  learning  augmented,  lateral-directional  flight  control  system.  Also  presented  in  this  report  is  the 
structure  of  the  computer  code,  instructions  for  running  the  program,  sample  input  and  sample  output. 
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TABLE  Bl.  BASIC  AIRCRAFT  PARAMETERS 


PARAMETER 

UNITS 

VALUE 

Gross  Weight 

lbs 

42097 

Mass 

slugs 

1308.4 

^xx 

12400 

Izz 

14300 

Ixz 

-2971 

S 

400 

b 

ft 

37.44 

FIGURE  Bl.  PERMISSIBLE  FLIGHT  ENVELOPE 


0  0.3  0.4  0.*  0.0  l.o  1.3  1.4  1.4  l.g  3.0 

Mach  »Mk«r 

PllMZSSXllZ 
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Pfdal  Porc«  -  lb 
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FIGURE  B2.  LATERAL  STICK  FORCE  VS.  STICK  DEFLECTION 

(REFERENCE  3) 


Ssj::  !a;or  Deflection 
deg'side 


Lateral  Stick  Deflection  -  in. 


Stabilator  Deflection  vs  Lateral  Stick  Deflection 


nCURE  B3.  PEDAL  FORCE  VS.  PEDAL  DEFLECTION 
(REFERENCE  3) 


DcflBCtton  •  in. 


DIRECTIONAL  PEDAL  FORCE  «•  PEDAL  DEFLECTION 

B-3 


DIFF  STABIUTOR  DEFLECTION  -  RADIANS 
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FIGURE  B4.  DIFFERENTIAL  STABILATOR  GEARING  CURVE 


RGUREBS.  AILERON  GEARING  CURVE 
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APPENDIX  C 

NONLINEAR  DATA  PLOTS 


C-1 


FIGURE  Cl.  ROLL  MODE  TIME  CONSTANT  (SEC) 
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FIGURE  C3.  DUTCH  ROLL  DAMPING 
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FIGURE  C5  VARIATION  OF  SIDE  FORCE  WITH  RUDOER  DEFL 

11/RAD) 
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FIGURE  C6.  VARIATIOIN  OF  ROLLING  MOMENT  WITH  RUDDER 

DEFL.  (1/RAD) 
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FIGURE  C7.  VARIATION  OF  YAWING  MOMENT  WITH  RUDDER 

DEFL.  (1/RAD) 
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FIGURE  CB.  VARIATION  OF  YAWING  MOMENT  WITH  YAW  RATE 

(1/RAD) 


FIGURE  eg.  VARIATION  OF  ROLLING  MOMENT  WITH  YAW  RATE 

(1/RAD) 
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FIGURE  Cll.  VARIATION  OF  ROLLING  MOMENT  WITH  AILERON 

DEFL.  (1/RAD) 
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FIGURE  C13.  VARIATION  OF  SIOEFORCE  WITH  SIDESLIP  (1/RAD) 
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FIGURE  CIS.  VARIATION  OF  ROLLING  MOMENT  WITH  SIDESLIP 

(1/RAD) 
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FIGURE  C16  VARIATION  OF  SIOE  FORCE  WITH  STABILIATOR 

DEFL.  (1/RAD) 
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FIGURE  C17.  VARIATION  OF  ROLLING  MOMENT  WITH  STABILATOR 

□EFL.  (1/RAD) 
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FIGURE  CIS.  VARIATION  OF  YAWING  MOMENT  WITH  5TABILIAT0R 

DEFL.  (1/RAD) 
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APPENDIX  D 
SAMPLE  INPUT  FILE 
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TABLE  Dl.  SAMPLE  FLIGHT  CONDITIONS 


CASE  NUMBER 

MACH 

ALTITUDE 

(FT) 

ANGLE  OF  ATTACK 
(RAD) 

1 

.2 

iiiiiii[iiiim]^iii[iiiiii 

.1414 

2 

.6 

35000 

.1093 

3 

.9 

10000 

.0171 

4 

0 

35000 

.05 

I^ote:  Case  4  is  not  a  valid  flight  condition. 


FIGURE  Dl.  INPUT  FILE  FORMAT 


Title 

Mach ,,A Ititude j(feet),Angle-Of-A ttack ;  (radians) 
Mach2,Altitude^(feet),Angle-Of-Attack2  (radians) 


Mach„Altitude„(feet),Angle-Of-Attack„  (radians) 
where  each  subscript  refers  to  each  flight  condition 


nGURE  D2.  SAMPLE  INPUT  FILE 

TEST  FILE  FOR  F-I8  TRANSFER  FUNCTION  PROGRAM  3/16/94 

.2, 1000, .1414 

.6,35000,. 1093 

.9,10000,.0171 

0,35000,.05 
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APPENDIX  E 
SAMPLE  OUTPUT  FILE 
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TEST  FILE  FOR  F-1 8  TRANSFER  FUNCTION  PROGRAM  3/16/94 
Case  Number  1 

mach  =  .200  alt  =  1000.0  alpha(rad)  =  .141 

Denominator  to  Lateral-Directional  Transfer  Functions. 

Coefficients  are  in  descending  order  of  s. 

1.000000  1.719688  3.037692  1.144893  .000000 

p  to  lateral  stick  force  transfer  function  numerator 
13.966980  17.499701  34.259322  .000000 

phi  to  lateral  stick  force  transfer  function  numerator 
13.966980  17.499701  34.259322 

beta  to  lateral  stick  force  transfer  function  numerator 

.000000  .000000 

r  to  lateral  stick  force  transfer  function  numerator 
.327676  .001629 

aycg  to  lateral  stick  force  transfer  function  numerator 
■7.843049  -8.566533  -19.231792 

p  to  pedal  force  transfer  function  numerator 
-.013788  -.000792  .111820  .000000 

phi  to  pedal  force  transfer  function  numerator 
-.013788 -.000792  .111820 

beta  to  pedal  force  transfer  function  numerator 
-.000385  -.068454 

r  to  pedal  force  transfer  function  numerator 
.068420  .000526 

aycg  to  pedal  force  transfer  function  numerator 
.006261  .000313  -.060769 

Case  Number  2 

mach  =  .600  alt  =  35000.0  alpha(rad)  =  .109 

Denominator  to  Lateral-Directional  Transfer  Functions. 

Coefficients  are  in  descending  order  of  s. 

1.000000  1.616000  2.561680  .967872  .000000 

p  to  lateral  stick  force  transfer  function  numerator 
29.709662  33.750175  59.906559  .000000 

phi  to  lateral  stick  force  transfer  fimction  numerator 
29.709662  33.750175  59.906559 
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beta  to  lateral  stick  force  transfer  function  numerator 

.000000 .000000 

r  to  lateral  stick  force  transfer  function  numerator 
.210279  .003787 

aycg  to  lateral  stick  force  transfer  function  numerator 
-16.683230-16.822641  -33.601716 

p  to  pedal  force  transfer  function  numerator 
-.001718  .000252  .306057  .000000 

phi  to  pedal  force  transfer  function  munerator 
-.001718.000252  .306057 

beta  to  pedal  force  transfer  function  numerator 
-.000113 -.072688 

r  to  pedal  force  transfer  function  numerator 
.072675  .000538 

aycg  to  pedal  force  transfer  function  numerator 
-.000177 -.000277 -.166413 

Case  Number  3 

mach  =  .900  alt  =  10000.0  alpha(rad)  =  .017 

Denominator  to  Lateral-Directional  Transfer  Functions. 
Coefficients  are  in  descending  order  of  s. 

1.000000  5.293200  18.167690  5.604424  .000000 

p  to  lateral  stick  force  transfer  function  numerator 
44.700194  221.409001  736.820090  .000000 

phi  to  lateral  stick  force  transfer  function  munerator 
44.700194  221.409001  736.820090 

beta  to  lateral  stick  force  transfer  function  numerator 
.000000  .000000 

r  to  lateral  stick  force  transfer  fimction  numerator 
2.462298  .047704 

aycg  to  lateral  stick  force  transfer  function  numerator 
-25.101046  -82.655994  -412.948207 

p  to  pedal  force  transfer  fimction  numerator 
-.050856  -.014647  3.596983  .000000 

phi  to  pedal  force  transfer  function  numerator 
-.050856  -.014647  3.596983 

beta  to  pedal  force  transfer  function  numerator 
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-.000668  -.266152 

r  to  pedal  force  transfer  function  numerator 
.265872  .001211 

aycg  to  pedal  force  transfer  function  numerator 
.017259  .003494-1.999366 

Case  Number  4 

mach  =  .000  alt  =  35000.0  alpha(rad)  =  .050 

Trim  flight  conditions  outside  of  permissible  flight  envelope 
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